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ABSTRACT 

Analytical procedures are developed for a composite system consisting of shape memory alloy fibers 
within an elastic matrix subject to uniform temperature fluxuations. Micromechanics for the calculation of the 
equivalent properties of the composite are presented by extending the multi-cell model to incorporate shape memory 
alloy fibers. A three phase concentric cylinder model is developed for the analysis of local stresses which includes 
the fiber, the matrix, and the surrounding homogenized composite. The solution addresses the complexities induced 
by the nonlinear dependence of the in-situ martensite fraction of the fibers to the local stresses and temperature, and 
the local stresses developed from interactions between the fibers and matrix during the martensitic and reverse 
phase transformations. Results are presented for a Nitinol/Epoxy composite. The applications illustrate the response 
of the composite in isothermal longitudinal loading and unloading, and in temperature induced actuation. The local 
stresses developed in the composite under various stages of the martensitic and reverse phase transformation are 
also shown. 


INTRODUCTION 

Shape memory alloys (SMAs) have found increasing applications as candidate smart materials and 
components of adaptive structures \ One potential application of shape memory alloys is their use as fibrous 
sensors or actuators within a composite material system. The development of such fibrous shape memory (SM) 
composites will provide weight savings, directionality in actuation and sensing, and some control over the resultant 
SM effect. Some technical details of such applications were discussed by Schetky and Wu 2 , Schetky 3 , and White 
et al 4 . The use of shape memory alloys as sensors may be related to the dependence of their properties on the 
martensite fraction (for example, their electric resistance) 5 . Conversely, the use of shape memory alloy (SMA) fibers 
as distributed micro-actuators is also very attractive because of the relatively high strains that can be generated. 

In summary, the shape memory (SM) effect is associated with a martensitic phase transformation of the 
material from a parent phase (ref. 1). During such martensitic transformations, the change in crystallography 
results in substantial strains, called transformation strains. Of particular interest, are alloys which can recover their 
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parent phase through a reverse phase transformation accompanied by recovery strains. These are SMAs which can 
undergo thermoelastic or complete phase transformations. In addition to the transformation strains, the phase 
transformations are accompanied by changes in the elastic, thermal and electric properties. The initiation of the 
martensitic or the reverse transformation depends on the free chemical energy in the material, which in turn, 
depends on the temperature and stress state. Hence, in a SM composite, the actuation process may be achieved by 
forcing the phase transformation to occur in the fibers by changing either the temperature or the stresses 3 in the 
fibers. 

Constitutive relations required to characterize the thermomechanical response of shape memory alloys 
were outlined in the works of Baumgart et al 6 , Tanaka 7 and Sato, Tanaka 8 and Brinson 9 . Applications of the 
unidirectional shape memory effect can be found in the papers of Liang and Rogers 10 , Lin and Rogers 11 , Baz et al 12 
and Baz and Chen 13 . Liang and Rogers 14 , Brinson and Lammering 15 and Boyd and Lagoudas 16 have provided 
formulations for the three-dimensional shape memory effect. 

The analysis of composites incorporating continuous shape memory alloy fibers seems a rather open-ended 
research area involving solution of several problems including admissible micromechanics and the nonlinear 
constitutive relations of SMAs. Indeed, the properties and the induced strain of the SMA fibers depend on the 
fraction of martensite, which in turn, depends on the in-situ temperature and effective stress state. Because of the 
fiber-matrix interaction in the composite, the developed in-situ stresses will depend on phase transformation strains, 
thermal strains and the constituent properties. This interdependence results in a nonlinear thermomechanical 
response (even by assuming thermoelastic austenite and martensite phases), which makes an iterative analysis 
seemingly necessary, although as shown in the present paper, an analytical solution can be obtained in some cases 
without iterations. Numerous micromechanical theories may be used to formulate the present problem, yet, limited 
work has been reported in the area of SM composites. Three-dimensional micromechanics were reported by Boyd 
and Lagoudas 16 using the Mori-Tanaka 17 micromechanics model. This model yields equivalent properties of 
composites that are in close agreement with results obtained from the finite element analysis, as was shown by Noor 
and Shah 18 . However, the assumption of average constituent stresses/strains may be an area of limitation, 
considering the dependence of martensite fraction on local stress and temperature fields. 

The present paper presents a combined micromechanical approach. A three-phase (SMA fiber, matrix, 
equivalent SM composite) concentric cylinder model is formulated for the accurate calculation of microstresses. The 
multi-cell micromechanics approach of Chamis and coworkers 19 is also extended to predict the properties of the 
composite medium encompassing the SMA fibers and the matrix. The solution is obtained for prescribed values of 
the effective stress in the fiber which yield the martensite volume fraction and the properties of the fiber and the 
composite medium, thus, avoiding an iterative procedure associated with readjustment of the properties at each step. 
Results quantify the isothermal loading-unloading response of the composite under various temperature and initial 
stress conditions, the capability of the composite to undergo temperature induced actuation, and the resultant local 
stresses in the composite during a loading-unloading cycle with stress induced phase transformation in the fibers. 


CONSTITUTIVE EQUATIONS 

The assumed constitutive equations of the SMA fibers, the matrix and the resultant SM composite are 
presented in this section. The fibers and the matrix are assumed to remain within the linear elastic range, which is a 
reasonable assumption since plasticity seems to usually occur after the martensitic transformation. Both the fibers 
and matrix are assumed isotropic. However, the composite is orthotropic with a plane of symmetry transversely to 
the fibers. 

Fibers- Assuming that both martensitic and parent phases are isotropic and that there is no directionality in 
the alignment of the martensite and austenite crystals, the SMA fibers may be assumed to remain isotropic, and their 
engineering properties are linear functions of the martensite fraction (rule of mixtures) 15 " 16 , 

(1) 
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where P represents the thermoelastic properties, subscript f indicates the fiber, £ is the current martensite volume 
fraction in the fiber, and superscripts M and A identify the martensite and austenite (parent) phases respectively. 

The behavior of a shape memory alloy fiber during the martensite or reverse phase transformation can be 
represented in a differential form 8 , 

<tf = [C f ] {e-cyf-co?} (2) 

where: overdot indicates time derivatives, s is the stress vector, e is the Green strain vector, T is the temperature, and 
x the martensite phase fraction; [C] and a are the current elasticity tensor and the vector of thermal expansion 
coefficients, respectively; co is the phase transformation tensor relating the transformation strains to the martensite 
volume fraction. The phase transformation tensor represents the recovery strain limits 10 . In the case where the 
tensors in eq. (2) are linear functions of x, eq. (2) can be integrated over time 10,15 to yield: 

o f -a° f = [C/Ol{e-a/OAr-Q?}- 


In thermoelastic martensitic transformations, the martensite fraction depends on the free chemical energy in the 
material, which in turn, includes the thermal and strain energy components. This relation may be represented by the 
following phenomenological equations 7,8 during the martensite transformation, 

§ = l-exp[^(7’"-7’) + c M aJ (4) 

and during the reverse transformation, 

5 r "- r K a d (5) 

T M S and T A S are the martensite and austenite phase start temperatures, o eff is the effective stress, the coefficients b M , 
c M , b A and c A are given by. 
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T M F and T a f are the martensite and austenite finish temperatures, and d M and d A are constants defining the 
slopes of the martensitic and reverse transformation lines, respectively. The characteristic temperatures depend also 
on the effective stress: 


3 



( 7 ) 


rr> M yi W nr A/ On A/ 

,/y. > - <ij ,if > 

<?/,?/> = <r/,r/> + 


+ 


J eff 


*M 


O 


eff 


where overbar indicates the start and finish temperatures of the respective phase transformation under reference 
stress-free conditions. Eqs. (4) and (5) may be rearranged to provide the range of the effective stress where stress 
induced martensitic and reverse phase transformations (at a prescribed temperature) will occur, respectively: 
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Note that in accordance with the metallurgical approach, eqs. (8) and (9) are obtained by identifying the 
values 5=1. and £=0.01 with the completion of the martensite and reverse transformations, respectively. 

Matrix . — The matrix is assumed to be isotropic and elastic, such that: 

= C m( € -«m* r ) (10) 


where AT is the difference between the current and stress-free temperatures, and subscript m indicates the matrix. 

Composite . — The equivalent constitutive equations for a composite medium with SMA fibers undergoing a 
martensitic or reverse transformation can be postulated to have the following form 

o c = [C c ]{e-a c r-to c ?} (11) 

where subscript c indicates the homogenized composite and x is the martensite phase fraction in the fibers. The 
main difference with eq. (2) is that the resultant composite will exhibit orthotropy in the LT plane (see Fig. 1) 
including the components of the phase transformation vector: 

= K’^r’VT (12) 

Within a range of linear dependence of the elastic, thermal and transformation components to the martensite volume 
fraction, integration of eq. (1 1) over time may yields, 

a c -a° c = C f (0[e-a e ©Ar-<o e q- 
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Both eqs. (3) and (13) remain valid during the martensitic and reverse transformation and they reflect the 
behavior of thermoelastic shape memory alloy materials stressed below their second yield limit associated with the 
plastic slip of the martensitic phase 20 . It is pointed out, that eq. (13) represents rather an approximation, because the 
transverse properties of the composite are not linear functions of x. In this case linearization over a range of x may 
be required, and eq. (13) should be applied within this range. In all cases, the initial conditions corresponding to the 
particular phase of the transformation should be considered, moreover, the continuity of stress-strain curves at the 
boundaries of adjacent transformation phases is preserved. It is finally emphasized that although equations (3) and 
(13) can describe both martensitic and reverse transformations, the corresponding stress-strain curves do not overlap 
each other, because of the differences between the effective start and finish stresses and temperatures defining the 
onset and completion of the martensite and reverse phase transformations [eqs. (4-5) and (8-9)]. 


MICROMECHANICS 


Multi-Cell Approach 


Good approximations of the equivalent properties of the previous SMA composite may be obtained using 
the multi-cell micromechanics approach 19 . According to the method, the representative cell is divided to a set of 
subregions (see Fig. la), each one containing portions of a square SMA fiber and the matrix. It is further assumed 
that the strains/stresses in the various regions of the representative cell are constant, the Poisson's effect on trans- 
verse deformations is negligible, and the fibers are packed in a square array formation. The analysis employs a 
mechanics of materials approach based on the equilibrium of stresses within each region and the compatibility of 
strains at the boundaries of each region. Thus, the equivalent incremental constitutive relations of the SMA 
composite are derived in accordance eq. (11). 

The resultant expressions for the elastic and thermal properties are similar to the ones reported for thermo- 
elastic composites, hence, are shown in the Appendix. According to the differential form of the constitutive relations 
for the fibers (2) and the composite (1 1), the phase transformation/recovery strain component has an analogous form 
to the thermal strain component. Therefore, at a given temperature and martensite fraction in the fibers, the 
equivalent transformation tensor of the composite w c may be evaluated similarly to the thermal expansion 
coefficient of the composite (see Appendix). 
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where, subscripts L and T indicate the longitudinal and transverse direction in the composite with respect to the 
fibers, k f is the fiber volume fraction, E f and E m are the moduli of elasticity of the fiber and matrix respectively, and 
E l and Ej. are the longitudinal and transverse moduli of elasticity of the composite medium. 


Concentric Cylinder Approach 

The previous multi-cell micromechanics is expected to yield sufficient approximations of the equivalent 
properties of the composite. Considering, however, the potential significance of local stresses (microstresses) in the 
response of the SMA fibers, a concentric cylinder model can provide better representations of the local stresses, 
which may have an effect on the axial response of the SMA composite. The three-cylinder assembly model shown 
in Fig. lb was formulated containing a fiber, surrounded by a matrix and an infinite SM composite medium with 
constitutive relations provided by either eqs. (1 1) or (13) and equivalent properties provided by the previously 
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described multi-cell micromechanics. The inclusion of the infinite composite cylinder is important, because it 
introduces the interactions between the fiber and matrix, and the surrounding composite due to the mismatch in the 
respective elastic, thermal and phase transformation components. Moreover, the consideration of a surrounding 
infinite composite, instead of a surrounding infinite matrix, extends the application of the approach beyond the low 
fiber volume ratio regime. 

Considering the axisymmetry of the problem, the equations of equilibrium within each material region of 
the concentric cylinder model (fiber, matrix, and composite) are reduced to the following equation: 


= o 

dr r 


(15) 


where subscripts x, r and 0 indicate longitudinal, radial and tangential directions, respectively. The axial strain is 
assumed to remain constant throughout the cell and the composite, while the tangential and radial strains are given 
by 



dw 
~d r 


(16) 


where w is the radial displacement. Eqs. (16) imply that the range of validity of the present solution is confined to 
geometrically linear deformations. 

Combinations of eqs. (15) and (16) and then substitution of constitutive eqs. (3), (10) and (13) for the fiber, 
matrix and composite regions respectively, yield an admissible solution in the form: 
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where A k j are constants of integration. The requirements of a bounded solution yield A, c = A/ =0. The remaining 
constants of integration are related via the conditions of continuity in the displacements and stresses at the fiber- 
matrix and matrix-composite interfaces: 


w f (r f ) = w m (r f ) 

drty = °7( r f ) 
» m (0 = » e (rj 


(18) 


where r f and r m are the radius of the fiber and the external matrix interface (Fig. lb), respectively. 

The effective stress, required for the determination of the martensitic fraction [eqs. (8-9)] and the state of 
the composite, is the following function of the fiber stresses: 

v eff = {o^-t/e) 2 +(o f r -o f d }} m (19) 


It is pointed out that the solution yields equal tangential and radial strains in the in the fiber, that is, 
e e =e=Aj f . Considering the isotropy in the transverse plane of the fiber, this also implies that the radial and 
tangential stresses in the fiber are equal and remain constant along the radial direction. This simplifies further the 
expression for the effective stress in the fibers: 
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(<V/ = l<Ml/ 


( 20 ) 


Computational Procedure 

As seen in the previous analytical formulation, the analysis of both the properties and the local stresses in a 
fibrous SM composite may require an incremental and iterative approach. This is due to the fact that the properties 
of the fibers, and, therefore, the average properties of the composite depend on the stresses via the martensite 
fraction. In turn, the stresses are affected by the properties of die constituent materials. To further complicate the 
matter, the start and finish temperatures of the martensitic and reverse transformation depend also on the stresses. 

However, the iterative process can be avoided and solution may be obtained in one step, within the range 
of x where eq. (13) is valid. In the later approach, the analysis begins with predefined values of the temperature and 
the effective stress in the fiber. Then the martensite fraction in the fibers is calculated according to eqs. (4) or (5). 
Once the martensitic function is defined, the properties of the fibers and the equivalent properties of the surrounding 
composite are calculated. Then, eqs. (15) and (20) are solved and the constants of integration A, f , A, m , A 2 m , A 2 C and 
the axial strain e x corresponding to the effective stress are found. The local strains and stresses are calculated. A 
final integration over the radius yields the average longitudinal stress in the composite 

/;- o x dr 


APPLICATION CASES 

The application cases presented in the following paragraphs illustrate the predicted stress-strain response of the 
composite under isothermal loading and unloading conditions, the capability of the composite to induce axial 
strains as a result of variations in the temperature, and the developed local stress fields in the material during stress 
induced martensitic and reverse transformations in the fiber. The shape memory material considered in the examples 
is nitinol Ni-44.8wt%Ti whose properties are presented in Table 1 as reported by Lei and Wu 21 . An epoxy matrix 
was assumed with properties shown in Table 1 and the fiber volume fraction was taken equal to 0.3. 


Composite Response 

The predicted loading and unloading cycle of the SM composite at 40°C is shown in Fig. 2. It is assumed that 
all local non-mechanical (shape memory and thermal) stresses are zero when the fibers are in the austenite phase 
(point A), that is, matrix consolidation occurs at 40°C with austenite fibers and the thermal residual stresses are 
zero. As seen in Fig. 2, the composite deforms linearly from point A to B because the fibers remain in the austenite 
phase (x=0). After point C the effective stress in the fiber becomes sufficiently high for the martensite transforma- 
tion to start and to induce transformation strains in the fibers, thus resulting in the nonlinear segment BC (0<x<l). 
At point C the fibers are fully transformed to the martensite phase (x=l) and the subsequent linear section CD 
illustrates the behavior of the composite with elastic martensite fibers. The different slopes of sections AB and CD 
represent the differences in the moduli of the martensite and austenite phases of the fiber. The predicted unloading 
phase is shown in section DEF. As the stress is reduced, the fibers remain in the martensitic phase until the reverse 
transformation starts (point E). Further stress reductions after point E, reduce the martensite fraction in the fibers 
and the recovery strains, resulting in the nonlinear curve EF (l>x>0). After point F the fibers are fully transformed 
to their parent phase (x=0) and the subsequent unloading follows the linear path FA. 
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The isothermal unloading of the previous SM composite system (reverse transformation in the SMA fibers) at 
four different temperatures is shown in Fig. 3. As in the previous case, the stress-free condition in the composite 
occurs with austenite fibers at 40°C. As seen in Fig. 3, at higher temperatures higher longitudinal stresses are 
required for starting and finishing the reverse transformation in the fibers. Fig. 3 implicitly illustrates the potential of 
SM composites to be used as strain or stress actuators by varying the temperature of the composite, as well as, it 
quantifies the ranges of stress, strain and temperature where such actuation is feasible. Within this range, the 
composite has the potential to generate significant longitudinal deformations or stresses. In all cases shown in Fig. 3, 
high recovery strains of an order of 6-8% have been calculated. Outside this range, the SM composite can not act as 
an actuator. However, as shown in the following case, it seems that this range may be modified by changing the 
processing conditions. 

Figure 4 also shows the isothermal unloading of the composite, however, matrix consolidation with fibers in the 
martensitic phase was assumed, such that the stress-free condition corresponds to fully transformed martensite fibers 
at 40 °C. In this manner, different states of local stresses are developed in the composite compared with the previous 
cases, that is, tensile longitudinal stresses in the fibers as the reverse transformation progresses. The predicted 
composite behavior differs significantly from the curves in Fig. 3. The reverse transformation occurs at different 
stress and strain ranges, and the recovery strains are smaller, of an order of 4% or less. This illustrates that the range 
of actuation of these novel composites may be modified by controlling the initial stresses in the composite. The 
predicted nonlinear unloading curve is also non-monotonic with respect to the stress. This is possible, because on 
the onset of reverse transformation tensile longitudinal and radial stresses will develop in the fiber (due to the 
negative transformation strain), which may suffice or exceed the values required to maintain a self-sustained reverse 
transformation. Hence, in order to achieve monotonic decrease in the martensite fraction (and the corresponding 
effective stress) in the fibers, the axial stress needs to increase and decrease again. This type of behavior is expected 
be observed in a "strain controlled" uniaxial test. 

The predicted capability of the composite to induce axial strain, under zero axial stress conditions, by 
controlling the temperature in the composite is shown in Fig. 5. The curve corresponds to a composite with stress- 
free austenite fibers. The curve begins at 20 °C when the fibers are in the martensitic phase. The reverse transforma- 
tion starts at 30 °C and finishes at 55 °C when the fibers have been fully transformed into austenite. It was predicted 
that, a 25 °C variation in temperature will result in a recovery strain of 6.3%. 


Local Stresses 

The predicted local stress distributions in the composite are presented here. The calculation of the stresses in the 
fiber and the matrix is important to ensure the integrity of the composite and prevent microcracking and fiber-matrix 
debonding. The longitudinal, radial and hoop stresses are shown in Figs. 6-9 as functions of the normalized radius 
for various phases of the isothermal stress-strain cycle shown in Fig. 2. The normalized radius in the figures is 
defined as the ratio of the local radius to the radius of the fiber (r/r f ). 

Figure 6 shows the stresses in the composite, when an average longitudinal stress s cx = -53.6 MPa is applied and 
the fibers still remain in the austenitic phase. Fig. 7 shows the stresses in the composite when the fibers undergo 
martensitic transformation (£ = 0.83) under the application of a tensile stress (s cx = 176 MPa). The applied tensile 
stress results in tensile longitudinal stresses in both the fiber and the matrix and in a tensile hoop stress in the matrix. 
Fig. 8 shows the local stress in the composite with fully transformed fibers ( x=l, s cx = 288 MPa). In spite the 
significant increase in the average applied stress, there is little change in the magnitudes of the radial and hoop 
stresses, but the longitudinal stresses in the fibers and the matrix increase significantly. 

Predicted stress distributions corresponding to the reverse phase transformation in the fibers ( x=0.40, s cx = 

12 MPa) are shown in Fig. 9. A comparison indicates that the stresses tend to return to a distribution similar to the 
one shown in Fig. 6. Further reduction of the longitudinal stress to -53.6 MPa resulted in full recovery of the parent 
phase in the fibers ( x=0) and yielded a stress distribution identical to the one shown in Fig. 6. This is not surprising, 
because at this point both loading and unloading curves coincide, hence, the corresponding internal stress states in 
the composite should be identical. 
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SUMMARY 


Micromechanics for the evaluation of the properties and analysis stress in unidirectional composite systems 
with continuous shape memory alloy fibers in an elastic matrix were presented. The multi-cell micromechanical 
method was extended to calculate the average properties of the composite when the fibers undergo martensite or 
reverse phase transformations. The axial response of the composite and the local stresses in the fiber and the matrix 
were determined using a concentric cylinder model, where the fiber and the surrounding matrix are embedded 
within an infinite composite medium. 

Evaluations of the mechanics were presented for a Nitinol/Epoxy composite system. The numerical results 
demonstrated the longitudinal response of the composite under isothermal loading and unloading conditions. They 
also illustrated the capability of the composite to induce substantial longitudinal strains under variation of temper- 
ature. The ranges of stress, strain and temperature where martensite and reverse phase transformations occur were 
also shown. It was predicted that substantial actuation strains may be induced in composites of intermediate fiber 
volume fraction (-0.3). It was also shown, that the range of stress and strain where the composite may be used as 
actuator may be modified by controlling the stress-free conditions in the composite. Finally, the distribution of 
microstresses was considered for an entire isothermal loading-unloading cycle of the composite material with stress 
induced phase transformation. The characteristic distribution of stresses remained unaffected, but their magnitude 
and sign changed drastically, depending on the martensite fraction of the fibers and the applied external load. In 
closing, this work has been the first step towards understanding and modelling the complex, yet promising, active 
response of continuous SMA fiber composites. Additional research will be conducted and reported in the future. 
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APPENDIX 


Elastic Properties of SM composites using the multi-cell method 19 : 


E L = E fr+ E m k m 


( 22 ) 



(23) 


G lt = G m 






1-Jkf 


/ G \ 

1 — =. 

G f) J 


(24) 


i-f, i~g 

(25) 

( /j 


= V A +V A 

(26) 

II 

to w 

? ■p 

t — » 

(27) 


Coefficients of Thermal Expansion: 

a L ~ + E m u m k m yE L (28) 



(29) 


where E m , G m , v m and a m are the normal and shear moduli of elasticity, the Poisson's ratio and the coefficient of 
thermal expansion of the matrix, respectively; and n f is the Poisson's ratio of the fiber which is independent of the 
martensite fraction. Volume fractions of the fiber and the matrix are denoted by k*- and k^, respectively. Subscripts L 
and T identify the longitudinal and transverse directions. 


11 



(Annealing temperature 600°C) 


Table 1 

Properties of Nitinol fibers 21 


Property 


Nitinol 

Epoxy 

E A 

(GPa) 

30.0 

3.4475 

e m 

(GPa) 

13.0 

• 

V 


0.33 

0.35 

a A 

(10- 6 /°C) 

12.5 

64.8 

a M 

(10- 6 /°C) 

18.5 

- 

1 S 

(°C) 

23 

- 

t m 

(°C) 

5 

- 

npA 

1 s 

(°C) 

29 


T A f 

(°C) 

51 


d M 

(MPa/°C) 

11.3 

- 

d A 

(MPa/'C) 

4.5 

- 

0) 


0.07 
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(a) Multi-Cell Model (b) Concentric Cylinder Model 

Fig. 1 Representative composite models, (a) Multi-cell; (b) Concentric cylinders. 
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Axial Strain 


Fig. 2 Isothermal stress-strain cycle of 30% FVR nitinol/epoxy composite (40 degs. C). Stress 
free conditions correspond to austenite fibers. 


20 degs C 
40 degs C 
60 degs C 
80 degs C 



Fig. 3 Isothermal unloading of 30% FVR nitinol/epoxy composite at various temperatures. 
Stress free conditions correspond to austenite fibers. 
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20 degs C 
40 degs C 
60 degs C 
80 degs C 


Fig. 4 Isothermal unloading of 30% FVR nitinol/epoxy composite at various temperatures. Stress-free 
conditions correspond to martensite fibers. 



Temperature, (degs C) 


Fig. 5 Temperature induced axial strain actuation with of 30% FVR nitinol/epoxy composite with 
reverse transformation in the fibers. Composite remains under stress-free conditions. 
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Fig. 6 Distribution of microstresses in nitinol/epoxy composite with fibers in the parent (austenite) 
phase ( ocx= -53.6 MPa , |=0). Isothermal loading cycle at 40 degs. C. 
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Longitudinal Microstress (Pa) 
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Fig. 7 Distribution of microstresses in nitinol/epoxy composite with fibers undergoing martensite phase 
transformation ( acx= 176 MPa, £=0.83). Isothermal loading cycle at 40 degs. C. 
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Fig. 8 Distribution of microstresses in nitinol/epoxy composite with fibers fully transformed to the martensite 
phase, (ocx=288MPa, £=1). Isothermal loading cycle at40degs. C. 
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